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Abstract 

To differentiate between the contribution of nnammary epithelial cells (MEC) and infiltrating innnnune cells to gene 
expression profiles of nnannnnary tissue during early stage nnastitis, we investigated in goats the in vivo transcriptional 
response of MEC to an experimental intra mannmary infection (IMI) with Staphylococcus aureus, using a non-invasive 
RNA sampling method from milk fat globules (MEG). Microarrays were used to record gene expression patterns during 
the first 24 hours post-infection (hpi). This approach was combined with laser capture microdissection of MEC from 
frozen slides of mammary tissue to analyze some relevant genes at 30 hpi. During the early stages post-inoculation, 
MEC play an important role in the recruitment and activation of inflammatory cells through the IL-8 signalling pathway 
and initiate a sharp induction of innate immune genes predominantly associated with the pro-inflammatory response. 
At 30 hpi, MEC express genes encoding different acute phase proteins, including SAA3, SERPINAl and PTX3 and factors, 
such as S100A12, that contribute directly to fighting the infection. No significant change in the expression of genes 
encoding caseins was observed until 24 hpi, thus validating our experimental model to study early stages of infection 
before the occurrence of tissue damage, since the milk synthesis function is still operative. This is to our knowledge the 
first report showing in vivo, in goats, how MEC orchestrate the innate immune response to an IMI challenge with 
S. aureus. Moreover, the non-invasive sampling method of mammary representative RNA from MEG provides a 
valuable tool to easily follow the dynamics of gene expression in MEC to search for sensitive biomarkers in milk for 
early detection of mastitis and therefore, to successfully improve the treatment and thus animal welfare. 



Introduction 

Mastitis is an inflammation of the mammary gland (MG) 
commonly caused by bacterial infection. Despite extensive 
management practices, it continues to be an economically 
important disease of dairy ruminants worldwide, due to 
reduced milk yield, milk discarded after treatment and 
cost of veterinary services [1]. A new challenge today is to 
reduce the use of antibiotics and treatments by increasing 
the natural ability of animals to resist infection. This strat- 
egy is critically dependent on a better understanding of 
the host immune response at the early stages of infection. 
Indeed, the establishment, persistence, and gravity of 
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infection are clearly dependent on the rapidity and effect- 
iveness of the host response against the invading patho- 
gen. Furthermore, clearance of bacterial pathogens from 
the gland is often governed by responses that occur within 
immediate hours after initial infection [2]. 

At the early stages of infection, the predominant de- 
fence strategy that is rapidly induced is the innate im- 
mune response. This response is ubiquitous, short acting 
and targets a range of different microorganisms [3-5]. 
Cells involved in this innate response are mainly repre- 
sented by infiltrating neutrophils and macrophages, 
which are the first cells to be recruited at the site of in- 
flammation [6]. They secrete cytokines, chemokines and 
use additional cellular defence strategies to kill invading 
bacteria [7,8]. However, little is known about the factors 
involved in this rapid recruitment at the site of infection. 
Hence, studies have increasingly focused on the mammary 
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epithelial cell (MEC) since it is the first cell to be con- 
fronted with the pathogen once it has entered the mam- 
mary gland. There are more and more lines of evidence 
indicating a prominent role for MEC in the initiation of 
the innate immune response that triggers subsequent neu- 
trophil infiltration into the mammary tissue, the activation 
of these innate immune cells, and an increase in milk 
somatic cell count [9]. 

However the MEC role has only recently been recog- 
nized essentially using in vitro approaches which allow 
examining the MEC specific contribution to the immune 
capacity of the udder without the immune cell expres- 
sion. The effect of different pathogens including Esche- 
richia coli, S, aureus [10], Streptococcus uberis [11], and 
bacterial cell wall components such as lipotechoic acid 
(LTA) or lipopolysaccharide (LPS) [12] were tested on 
different mammary cell lines in either primary isolates of 
bovine MEC (pbMEC) or MAC-T cells, an SV40 immor- 
talized bovine MEC line [12,13]. These studies indicate 
that MEC respond robustly and rapidly to challenges 
with low levels of bacteria or bacterial cell component 
and show that MEC is able to generate a variety of in- 
flammatory mediators such as cytokines, chemokines 
and host defence peptides (p-defensins). However, these 
reports underlined some limits associated with in vitro 
experiments and show differences between mammary 
cell cultures (pbMEC and MAC-T) [12]. Furthermore, 
cells are studied out of their physiological context and 
consequently, they do not properly reflect changes in 
gene expression induced by mastitis in the udder. Thus, 
a disproportionate increase in chemokines such as CCL5 
has been shown between pbMEC and the whole mam- 
mary gland [14]. In addition, there is only one study pro- 
filing the extents of global changes in the transcriptome 
of pbMEC during the early stages of infections with 
E, coli and S, aureus, reflecting in part the kinetics as- 
pects of immune response regulation [15]. Consequently, 
the role of MEC in the initiation of the innate immune 
response, in their physiological context and in the very 
early stages of infection, remains poorly defined. 

Recently, we have reported [16] that substantial quan- 
tities of high-quality RNA, coming from MEC, can be iso- 
lated from milk fat globules (MFG). This non invasive 
method allows easy and repetitive sampling without dam- 
aging mammary tissue, providing a significant improve- 
ment and a valuable tool to assess gene expression of the 
mammary secretory epithelium during the course of an in- 
fection. Furthermore, with the availability of microarray 
technology, which makes it possible to examine complex 
interactions between host and bacterial pathogens [17], it 
has become feasible to determine the expression of mul- 
tiple genes simultaneously. We successfully applied this 
approach to follow the response of MEC early in the infec- 
tion process of the goat mammary gland. 



We report here that MEC rapidly react (a few hours 
after the experimental infection) by secreting the chemo- 
tactic factor interleukin-8 (IL-8), which is known to be 
one of the major mediators of the inflammatory re- 
sponse involved in the recruitment and activation of 
immune cells, including neutrophils. Furthermore, we 
provide strong evidence that MEC specifically expressed 
several acute phase proteins (APP) such as serum amyl- 
oid A3 (SAA3), pentraxin 3 (PTX3) and alpha- 1 antipro- 
teinase (SERPINAl). This work also indicates that the 
production of different APP by MEC, easily detectable 
from milk in the early stages of IMI, could provide sensi- 
tive biomarkers for early detection of mastitis and there- 
fore, to successfully improve its treatment and thus 
animal welfare. 

Materials and methods 

Animals 

Five healthy Alpine goats, at the early peak of lactation 
(30-40 days of milking), at first parity and without intra 
mammary infection (IMI) were selected and housed at 
Centro Zootecnico Didattico Sperimentale of Facolta of 
MedicinaVeterinaria (Universita degli studi di Milano). 
Goats were monitored for IMI throughout the lactation 
period with 3 weekly half udder milk sample analysis be- 
fore challenge. Bacteriological analyses were based on 
procedures previously described [18]. Briefly, ten microli- 
ters of each milk sample was spread on blood agar plates 
(5% defibrinated sheep blood). The plates were incubated 
aerobically at 37 °C and examined after 24 h and 48 h. 
The colonies were provisionally identified based on Gram 
stain, morphology, and hemolysis pattern, and the num- 
bers of each colony type were recorded. The representa- 
tive colonies were then subcultured on blood agar plates 
and incubated aerobically at 37 °C for 24 h to obtain pure 
cultures. Catalase and coagulase production was tested for 
gram-positive cocci. Specific identification of staphylo- 
cocci was made using commercial micro-methods (API 
Staph; BioMerieux, Italy). The infection status of milk 
samples was defined according to the procedures recom- 
mended by the National Mastitis Council (NMC, 1987) 
and IMI was diagnosed when > 500 cfu/mL and 1 to 3 col- 
ony types were isolated. Milk samples from which many 
colony types or < 500 cfu/mL of any microorganism 
were isolated were regarded as contaminated or unin- 
fected, respectively. At the moment of challenge, no in- 
fection in the udders was observed, as was confirmed by 
the absence of mastitis pathogens in foremilk samples 
tested for three consecutive days just before the experi- 
mental challenge. 

Experimentally induced mastitis and milk samples 

Prior to intra mammary challenge, both half-udders 
were milked by hand and the teat ends were carefully 
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disinfected with chlorhexidine. The left half-udder of each 
goat was inftised with 1 mL (10^ CFU/mL) inoculum of 
Staphylococcus aureus strain DV137 which was originally 
isolated from a chronic case of caprine mastitis. The right 
half-udder, used as an uninfected control, was infused with 
1 mL of sterile pyrogen-free PBS (Figure 1). 

First, milk samples from uninfected right (Rx, TO) and 
left half-udder (Lx, TO) were collected before IMI with 
S, aureus. Three hours later, intra mammary challenge 
with inoculum of S, aureus was carried out. Subsequent 
milk sample collections were carried out at 6, 12, 18, 24 
and 30 h after challenge from both right (uninfected) 
and left (infected) half-udder, taking care to avoid ribo- 
nuclease (RNase) contamination: cleaning the udder and 
teats with a clean linen impregnated with an antiseptic 
solution followed by a spray of RNAse Zap, drying with 
a paper towel (Kimwipes) and performing a manual 
milking with disposable nitrile gloves. For each animal 
and time point, 150 mL of milk were collected from 
each half-udder into sterile, RNase free tubes (3 x 50 mL 
Falcon tubes). Samples were immediately kept on ice 
prior to MFG collection. 

The clinical scoring system to classify mastitis symp- 
toms cases was as described for dairy cows [19]. In 
addition to abnormal milk, this system is based on 
measurement of rectal temperature, hydration status and 
clinical attitude. Severity of clinical signs was scored as 
mild, moderate or severe. A mild score was assigned 
when the milk was grossly abnormal and no other local 
or systemic signs of inflammatory disease were seen; a 
moderate score was assigned when the milk was grossly 
abnormal and there was firmness or swelling of the af- 
fected mammary gland, but none or only one of the sys- 
temic signs of inflammatory disease are seen. A severe 



score was assigned if the milk was grossly abnormal, 
there was firmness or swelling of the affected mammary 
gland, and at least 2 of the following systemic disease 
signs were seen: rectal temperature > 39.5 °C, hydration 
score showing moderate to marked enophthalmos, and 
attitude score showing signs of marked depression [19]. 

All experimental procedures were performed accord- 
ing to the Italian legislation, following approval by the 
ethics committee of University of Milan. 

S. aureus bacterial counts 

For determination of S, aureus bacterial counts (cfu/mL), 
series of dilutions (10^ to 10^) were prepared from 1 mL 
of each sample of milk, diluted with 9 mL of a 0.1% saline 
peptone solution. Then 0.1 mL of the serial dilutions were 
inoculated on the surface of Baird-Parker agar and spread 
with a spatula. The incubation was done at a temperature 
of 37 °C for 16 h (overnight incubation). 

Milk fat globule collection 

Samples were centrifuged at 2000 x g for 10 min at 4 °C 
to isolate milk fat. The supernatant fat layer was trans- 
ferred to a new 50 mL-Falcon tube using a sterile spatula. 
Then, 500 (iL of fat were put into a 15 mL-Falcon and 
1.5 mL of TRIzol® LS solution (Invitrogen, Life Technolo- 
gies, Carlsbad, California, USA) was added and the tube 
was vortexed vigorously prior to storage at -80 °C. The 
entire process of milk sample collection and storage of 
MFG was completed within 2 h and all procedures were 
carried out at 4 °C. 

Tissue collection for laser capture microdissection 

At the end of the experimental protocol (30 hours post 
infection (hpi)), goats were slaughtered, according to 
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Figure 1 Experimental infection worlcflow. The left udder was challenged by 5. aureus whereas the right udder remained uninfected as the 
control. Then every six hours, milk was sampled and centrifuged to extract RNA from Milk Fat Globule (MFG). At 30 hpi, goats (n = 5) were 
slaughtered and mammary tissue samples were taken for MEC capture experiments using LCM (total RNA extraction from micro-dissected MEC) 
and total RNA extraction from deep alveolar mammary parenchyma (referred to as MG). Finally, these different sources of RNA were analyzed 
using microarrays (only MFG) and qPCR. 
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surgical and experimental procedures in compliance with 
the policy of INRAs Animal Care Committee, after milk 
sampling. Tissue samples were collected aseptically ft-om 
the five goats within 10 min after slaughtering. A piece of 
deep alveolar parenchyma, without visible connective tis- 
sue, was removed from the left udder (infected) and right 
udder (uninfected). The collected tissue was washed in 
cold PBS solution (on ice), 5 mm^ pieces of tissue were cut 
and embedded into OCT® (TissueTek™) in a cryomold of 
1 cm^ (Bayer™) and immediately placed on dry ice or in a 
SnapFrost™ system (Alphelys, Elancourt, France) contain- 
ing isopentane at -80 °C. Samples were stored at -80 °C 
until ftirther processing. The time delay between slaughter- 
ing and tissue freezing was less than 20 min. Laser Capture 
Microdissection (LCM) was carried out using the Veritas 
Microdissection system and software (Arcturus, Life Tech- 
nologies, St Aubin, France), as previously described [20]. 

Total RNA extraction 

Total RNA was extracted from mammary tissue samples 
(pieces of deep alveolar parenchyma referred to as MG) 
taken at 30 hpi on the left (infected) and the right (unin- 
fected) half-udder, using TRIzol Reagent (Invitrogen, Life 
Technologies) according to the manufacturers instruc- 
tions. Total RNA was extracted from MFC using TRIzol® 
LS solution (Invitrogen, Life Technologies) following the 
original manufacturers protocol with slight modifica- 
tions, essentially as described by Brenaut et al. [16]. Prior 
to RNA quality control, a DNase treatment was carried 
out to remove any contaminating genomic DNA accord- 
ing to the manufacturer s protocol (Qiagen, Courtaboeuf, 
France). 

Total RNA was extracted from captured cells using the 
PicoPure® RNA Isolation Kit (Arcturus, Life Technologies) 
according to the manufacturers instruction protocol, 
including on-column RNase-free DNase I treatment 
(Qiagen). CapSure macrocaps with captured cells were 
inserted into RNase-free 500 [iL microcentrifuge tubes 
containing 25 [iL of extraction buffer (XB). The tubes were 
inverted to allow the reaction between the buffer and the 
surface of the cap. RNA were extracted from scraped sec- 
tions (tissue remaining on the slide after capture) by pipet- 
ting 50 [iL of XB buffer onto the remaining tissue on the 
glass slide, which was then gently scraped off and trans- 
ferred into RNase-free 500 [iL microcentrifuge tubes. 
RNA from cap and section scrapes were eluted with 15 (iL 
and 30 [iL of elution buffer (EB), respectively. 

RNA quality control and single strand cDNA synthesis 

Purity, concentration and integrity of total RNA intended 
for microarray analysis and qPCR were assessed using 
two independent techniques. RNA purity was evaluated 
by absorbance readings (ratios A260/A230 and A260/ 
A280) using the NanoDrop ND-1000 spectrophotometer 



(Thermo Fisher Scientific, Villebon-sur-Yvette, France). 
The fluorimetric method and micro-capillary electrophor- 
esis device developed by Agilent Technologies (Les Ulis, 
France) was used to determine RNA concentration and 
quality with RNA 6000 Pico LabChip" Kit in the Agilent 
Bioanalyzer 2100 system. Quality was evaluated using 
the RNA Integrity Number (RIN) value, introduced by 
Schroeder et al. [21]. A PCR test was routinely performed 
to ensure that RNA samples were free from genomic 
DNA (primers matching with intron sequences). First- 
strand cDNA was synthesized from 5-10 ng total RNA 
primed with oligo(dT)2o and random primers (3:1, v/v) 
using Superscript III reverse transcriptase (Invitrogen, Life 
Technologies) according to the manufacturers instruc- 
tions. Then, 1 [iL RNase H (2 U/(iL, Invitrogen, Life Tech- 
nologies) was added to the reaction mix which was 
incubated for 20 min at 37 °C to remove RNA from het- 
eroduplex. cDNA, thus obtained, was stored at -20 °C. 

Microarray processing and data analysis 

A microarray analysis was performed, in a mono color 
experimental design, to compare gene expression pro- 
files of MFC at 4 time points of infection (0, 12 h, 18 h 
and 24 h) on the 5 goats. RNA samples were labelled 
using T7 RNA polymerase, which simultaneously ampli- 
fies target material and incorporates cyanine 3-labelled 
CTP (Agilent Technologies) according to the manufac- 
turers instructions. Briefly, total RNA (120 ng) was first 
reverse transcribed and then converted to labeled cRNA 
with Cy3 dyes. Before hybridization, quantity and quality 
of purified cRNA was assessed with the Agilent 2100 
Bioanalyzer. 

Since goat microarrays were not available, we used the 
Sheep Gene Expression Microarray, 8 x 15 K (Agilent 
Technologies) developed for a closely related species 
{Ovis aries), to profile gene expression. Hybridization 
was performed according to the manufacturer s instruc- 
tions. Briefly, each microarray was hybridized with 
fluorescent labelled (Cy3) cRNA samples. Three hun- 
dred ng of each labelled cRNA were fragmented for 
30 min at 60 °C in 25X fragmentation buffer. Then, 
fragmented and labelled cRNA samples were diluted in 
2X GE X Hybridization Buffer Hi-RPM and hybridized 
onto the 15 K ovine microarray for 17 h at 60 °C in a ro- 
tating hybridization oven. After hybridization, microar- 
rays were washed twice with washing buffer 1 at room 
temperature for 1 min and then with washing buffer 2 at 
37 °C for 1 min. Microarrays were scanned on an Agilent 
scanner G2565BA using the extended dynamic range scan 
mode. Data were processed (including Lowess normalisa- 
tion) and extracted with Feature Extraction software ver- 
sion 10.5 (Agilent, Technologies). 

Statistical analysis of microarray data was carried out 
using the commercially available software package 
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GeneSpring® GX vlO.0.2 software (Agilent Technolo- 
gies), using the latest gene annotations available on a 
web-accessible resource established by SIGENAE [22], 
Information about this experiment has been deposited 
in NCBI Gene Expression Omnibus (GEO) with accession 
number (GSE39315). All data are "Minimum Information 
About a Microarray Experiment" (MIAME) compliant. 
First, we applied a quantile normalization on data, which 
assumes that the distribution of gene abundance is the 
same for all samples, and spots that did not meet minimum 
signal intensity were removed. The resulting signal infor- 
mation was analyzed using repeated measures ANOVA 
with a Benjamini-Hochberg correction for multiple com- 
parisons. This approach was used since the data correspond 
to a longitudinal study which measures each individual at 
each of several time points. We note that a standard 
ANOVA assumes independence between time points and 
is not able to differentiate between the variability within 
each individual (i.e. across time points) and the variability 
between individuals, resulting in inflated error variance 
estimates. By accounting for the correlation among re- 
peated measurements for each individual within the 
model, these error variance estimates can thus be 
reduced. This analysis included pair wise comparisons 
between time points (0 h vs. 12 h, 12 h vs. 18 h, and 
18 h vs. 24 h), which yielded a list of probes that were 
differentially expressed in at least one of the pairs con- 
sidered. A first filter was applied to select the genes that 
displayed an adjusted j^-value less than 0.05. The output 
of this analysis was additionally filtered by fold expres- 
sion, generating lists of differentially expressed genes 
with at least 1.5-fold change. This final set of genes was 
then referenced to biological functions with Ingenuity 
Pathways Analysis software v6.0 (IPA, Ingenuity Systems, 
Redwood City, CA, USA). The contribution of genes to 
the identified function was explained with an associated 
P'Vdlue, calculated with a right-tailed Fisher s exact test. 

To identif)^ temporal relationships among this subset 
of genes, a gene regulatory network was reconstructed 
from the expression data using an approach specifically 
developed for replicated longitudinal data, called Empir- 
ical Bayes Dynamic Bayesian Networks (EBDBN), imple- 
mented in the ebdbNet package in R [23]. Briefly, the 
EBDBN method makes use of an iterative empirical 
Bayes estimation procedure for a linear Gaussian state- 
space model. This approach has the advantage of being 
able to simultaneously model linear relationships among 
genes and a set of unobserved hidden variables (e.g., 
proteins, transcription factors, or other unobserved cel- 
lular entities) over time, as well as to explicitly model 
feedback loops of expression from one time point to an- 
other. Using the approach detailed in the package, the 
dimension of the hidden variables was chosen to be 2. 
The model was run for 10 different initializations, and 



gene-to-gene interactions were retained only if identified 
in at least 9 runs at a significance level of 95%. 

Primer design and real time quantitative PGR 

RNA isolated from goat MFG and microdissected MEC 
were used to determine the expression of selected genes 
of interest by real time quantitative PGR (qPCR) as de- 
scribed by Bevilacqua et al. [24], using the SYBR green 
PGR master mix and an ABI Prism 7900 Sequence 
Detection System (Applied Biosystems, Life Technolo- 
gies). The reference sequence for primers was taken 
from a bovine database. Expressed sequence tags (EST) 
of caprine sequences found by the BLASTn search and 
reference were aligned to create one consensus. To elim- 
inate the risk of genomic DNA amplification we system- 
atically chose primers that hybridize on exon-exon 
junctions. Only junctions with great homology between 
sequences were chosen for primer design. Then primers 
were designed using Primer Express Software version 2.0 
(Applied Biosystems, Life Technologies), and purchased 
from MWG Biotech (Table 1). First, primer efficiency 
was validated with a standard curve of four serial dilu- 
tion points (ranging between 1 ng and 1 pg of reverse 
transcribed total RNA) and a no template control 
(NTC). A qPCR amplification mixture (20 [iL) contained 
5 [iL single strand cDNA template, 10 [iL 2X Power SYBR 
Green PGR Master Mix buffer (Applied Biosystems, Life 
Technologies) and 1.2 [iL of forward and reverse primers 
(5 [iM) to reach a final primer concentration of 300 nM. 
After optimization of qPCR systems (efficiency ranging 
between -3.30 and -3.45), amplification reactions were 
run (in triplicate). Tests for non-amplification of genomic 
DNA was carried out systematically. The results generated 
by the Sequence Detection Software (Applied Biosystems, 
Life Technologies, version 2.3) were exported as tab- 
delimited text files and imported into Microsoft Excel for 
further analyses. Relative quantification analysis was per- 
formed using the software program qBase [25] in which 
relative expression levels were normalized with respect to 
the selected reference genes {RPS24 and PPIA), following 
the "Minimum Information for Publication of Quantitative 
Real-Time PGR Experiments" (MIQE) guidelines [26] . 

Statistical analysis 

Reliability of reference genes {RPS24 and PPIA) was 
evaluated with the GeNorm Visual Basic application for 
Microsoft Excel as described by Vandesompele et al. 
[27]. Statistical analyses were performed using the R 
software, with the Rcmdr package. The non-parametric 
Mann- Whitney U-test was used to analyse the differ- 
ences in relative expression for each gene of interest be- 
tween MG and MFG, non infected versus infected, at 
12 h, 18 h and 24 hpi but also between infected MG and 
microdissected MEC, at 30 hpi. 
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Table 1 Primer sequences used for qPCR experiments 



Genes 


Primers 


Sequence 5' > 3' 


Amplicon size, nt 


as2-casein iCSNlS2) 


forward 


CTG G^ ATG G^ GGA CTG GAA AA 


76 




reverse 


AAC ATG CTG G^ GTA TGA AGT AAA GTG 




K-casein {CSN3) 


forward 


AGG TGC AAT GAT GAA GAG m m C 


66 




reverse 


CCC AAA AAT GGC AGG G^ AA 




Cluster of differentiation 3 epsilon (CD3e) 


forward 


ACG CTG TAC CTG AAA GCA AGA 


118 




reverse 


AAT ACA CCA GCA GCA GCA AG 




Cluster of differentiation 68 {CD68) 


forward 


GAT CTG CTC TCC CTG AAG CTA CA 


79 




reverse 


CAT TGG GAC AAG AGA AAC TTG GT 




Cluster of differentiation 18 (CD/8) 


forward 


AGC GAC CTC AGG GAG TAC CA 


65 




reverse 


™ TCG TTG JTC CAC TGG GAC TT 




Chemokine 4 (CC/.4) 


forward 


CAG CCG TGG TAT TCC AGA CC 


109 




reverse 


CTC GGA GCA GCT CAG TTC AGT 




Cyclophilin A {PPAI) 


forward 


TGA C^ CAC ACG CCA TAA TGG T 


62 




reverse 


CAT CAT CAA A^ TCT CGC CAT AGA 




Galectin 3 


forward 


GTG GTA AAC C^ TCA AAA TAC AAG TGC 


101 




reverse 


A^ TTT CAC CCG ATG A^ GTA CTG 




G protein-coupled receptor 97 {GPR97) 


forward 


GAG ATC ACC JTC TCC CAC CAG 


204 




reverse 


TGT GGA GCA GCC CAA GGA 




Interleukin 8 ilL-8) 


forward 


TGA GAG TGG GCC ACA CTG C 


103 




reverse 


CAC AAC C^ CTG CAC CCA C^ 




Keratin 14 {KRT14) 


forward 


CCC AGC TCA GCA TGA AAG C 


57 




reverse 


AGC GGC C^ TGG TCT C^ C 




Serum amyloid A3 iSAA3) 


forward 


CTG GGC TGC TAA AGT GAT CAG TAA C 


69 




reverse 


CCC TTG AGC AGA GGG TCT GTG A^ 




S100A12 


forward 


TCC ACC AGT ACT CCA TCC GG 


102 




reverse 


TGG TGT TTT TGA GGC AGT TGG 




Serglycin 


forward 


TCC AGC AGA ATC CCA CCT CTA 


107 




reverse 


CCA GAA CCT GAT CCT GAG ACG 




SerpinAl 


forward 


AAG AAA TAT GCA AGT TCT GCC AAT T 


101 




reverse 


ACC CTG TTG ATG CCC AGT TC 




TNFa 


forward 


CAG AGG GAA GAG CAG TCC CC 


101 




reverse 


TGG GCT ACC GGC TTG TTA TTT 




Toll-like receptor 2 aLR2) 


forward 


TAA ACT TGA GAG TGG AGG TCA AAT CA 


101 




reverse 


TCA GAG GCT CCT TCC GTG G 




iiiLeiieuKin lueia (^iL-ipj 


fo rwa rd 


r-,Ar AAT AAr; att cct mr, err 


1 ni 

1 U 1 




reverse 


TCT ACT TCC TCC AGA TGA AGT GT 




Pentraxin 3 (PTX3) 


forward 


CCG AGC TGT GCA GGG CT 


101 




reverse 


GCA CGC TTG CAA AAA TCT TCT T 




Cathelicidin 


forward 


GAG AAT GGG CTG GTG AAA CAG 


107 




reverse 


GGG CGA AGT CTC CTC ACA CTC 




24S ribosomal protein {RPS24) 


forward 


TGG TGGTGG CAA GAC AAC TG 


66 




reverse 


TTC TTC GCG TAA TCC AAG GAA 





Each pair of primers amplifies the target cDNA (amplicon sizes ranging between 57 and 248 nucleotides) In Its 3' UTR. Primer pairs were designed with primer 
Express Software v2.0 (Applied Biosystems) except for 24S ribosomal protein primers, which were manually designed. 
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Results 

Kinetics of pathological features: an increase of S. aureus 
in milk followed by an increase of somatic cells in milk 

Prior to experimental infection, goats were controlled to 
be free from S. aureus and other IMI throughout the 
three weeks of lactation to avoid the influences of adap- 
tative immune mechanisms. Following administration of 
10^ CFU/mL of S. aureus DV137 to the left udder and 
the same volume of sterile pyrogen-free PBS into the 
right udder, left udders contained only the pathogen 
used for infection whereas the control udder (right) 
remained bacteriologically negative during the observa- 
tion period. A high number of pathogens were found 
18 h after infection, reaching 4 x 10^ CFU/mL and then 
decreasing to 1.5 x 10^ CFU/mL at 24 hpi (Figure 2). 
Mild mastitis symptoms were concomitant with in- 
creased somatic cell counts (SCC) in milk. In infected 
udders, the number of somatic cells significantly increased 
from 893.8 x 10^ at 18 hpi to 53 x 10^ cells/mL at 24 hpi, 
while SCC in PBS -infused udder halves showed no signifi- 
cant change throughout the period with a mean value of 
980 (± 147) X 10^ cells/mL (Figure 2). Thus, all infected 
goats showed a sharp increase of S, aureus in milk at 18 
hpi which was followed by a significant increase of SCC in 



Kinetics of infection 




B 

Kinetics of somatic cell recruitment 

12000000 1 




0 6h 12h 18h 24h 



Figure 2 Kinetics of infection (A) and somatic cell recruitment 
(B). A: Staphylococcus aureus colony forming unit (CFU/mL) 
evolution during the 24 first hpi. B: Somatic cell counts (SCC/mL) in 
the left (infected) udder of the 5 goats challenged by 5. aureus. The 
control PBS-infused udders remained free from detectable infection 
throughout the study for the five goats. 

v J 



milk at 24 hpi. Thirty hours after experimental challenge 
all goats showed severe signs of clinical mastitis: the milk 
was grossly abnormal, there was firmness and swelling 
of the affected left mammary gland and the rectal 
temperature was > 39.5 °C (40.2 ± 0.5 °C). 

RNA sampling from MEC in their physiological context, 
during infection 

Two different sampling methods were used to access 
MEC materials from each of the 5 goats throughout the 
time course of the IMI challenge with S, aureus (Figure 1). 
The first method consisted in milk collection, every six 
hours post infection, and extracting RNA from cytoplas- 
mic crescents trapped in MFC during the secretion 
process from MEC. The RNA integrity number (RIN) of 
total RNA isolated from MFG was estimated, on average, 
to 7.5 ± 0.5. Then, at the end of this first part of the proto- 
col, mammary tissues were collected on slaughtered goats 
{n = 5) and MEC captured from sections of mammary tis- 
sue by Laser Capture Microdissection (LCM). RIN for 
RNA extracted from microdissected cells were of compar- 
able quality (6.8 ±0.4) with RNA obtained from MFG. 
Combining these two techniques is a convenient approach 
to follow MEC response in a physiological context, pro- 
vided that contamination by other cell types can be 
excluded. To assess such an eventuality, qPCR systems 
designed to quantify specific markers for mammary myoe- 
pithelial cells {Krtl4), lymphocytes {CD3e), neutrophils 
{CD18), and macrophages {CD68 and CD18) were imple- 
mented [16]. 

RNA extracted from MFG highly reflect MEC gene 
expression pattern, up to 18 hpi 

During MFG preparation, some immune cells (mainly 
neutrophils and to a lower extent macrophages) present 
in milk can be trapped in the cream layer. To evaluate 
the level of RNA contamination, qPCR analyses were 
carried out using specific gene markers of these immune 
cells: CD68 for macrophages and CD18 for neutrophils 
and macrophages. Two markers of MEC {CSN3, CSN1S2) 
were also included as positive controls which were 
expressed at the same level in MFG of uninfected and in- 
fected udder, at all time points post infection (Figure 3). 
Regarding the level of contamination by macrophages and 
neutrophils, qPCR analyses showed that these immune 
cells were actually present but in very low amounts, be- 
tween 7 to 9-folds lower in MFG than in the whole 
deep mammary alveolar parenchyma (MG), until 18 hpi 
(Figure 3). These very low values (RQ mean CD68 = 
0.15 ±0.05 and RQ mean CDi8 = 0.17 ± 0.03, between 
TO and T18) reflect a slight contamination by macro- 
phages and neutrophils, during MFG preparation. Con- 
tamination increased and became non negligible, 
essentially for macrophages, with an RQ mean for 
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Figure 3 Determination of specific markers by qPCR to assess the level of contamination of RNA extracted from milk fat globules by 
RNA from immune cells. Relative expression (± SEM) is given for specific gene marl<ers for MEC CSN3 (K-casein, green), CSN1S2 (as2-casein, blue), 
macrophages {CD68, red) and macrophages + neutrophils (CD/8, orange) in the left half-udder (infected) or the right half-udder (uninfected), 
sampled at 30 hpi on slaughtered goats (mammary gland, left), and in milk fat globules (right), before infection (uninfected) and after IMI challenge 
with S. aureus (12 hpi, 18 hpi and 24 hpi). ""Significance is relative to fold change in expression of uninfected MG (adjusted p-value < 0.05). 



CD68 reaching 4.96 ± 0.31, at 24 hpi, whereas it was 
30-folds lower at 18 hpi. Therefore, at least until 18 
hpi, responses to IMI challenge observed via MFG tran- 
scriptome analysis have to be essentially attributed to 
MEC, whereas from 24 hpi and beyond, a part of the 
gene expression data is influenced by infiltrating im- 
mune cells trapped in the cream layer during MFG 
preparation. 

Selectivity of MEC capture by LCM on cryo-sections of 
mammary tissue 

To provide additional insight into the immune role asso- 
ciated with MEC during infection, we investigated gene 
expression profiles of microdissected MEC. Indeed, 
24 hpi, MFG extracted material was no longer represen- 
tative of MEC owing to a significant contamination of 
MFG by immune cells (see above). Therefore, at the end 
of the experimental milking sampling protocol (i.e. 30 
hpi), mammary tissue was collected on slaughtered goats 
in order to capture MEC from tissue cryo-sections, 
using LCM (Figure 1). This technique makes it possible 
to capture MEC while they are in their physiological 
context and therefore to obtain RNA specifically from 
this cell type at 30 hpi, thus allowing a comparison with 
RNA isolated from MFG at the same stage. However, in- 
fection leads to an infiltration of immune cells within 
the secretory parenchyma which affects tissue morph- 
ology and reduces the surface on which LCM could be 
carried out (Figure 4B). Consequently, only four goats 
were sampled given the morphology which was strongly 



affected by infection and therefore not easily suitable for 
microdissection. To assess the selectivity of MEC cap- 
ture, mRNA transcript level of a MEC specific marker 
{CSN3) was estimated by measuring the relative expres- 
sion between captured cells and their corresponding 
mammary tissue scrapes. The amount of mRNA encoding 
CSN3 was unchanged (RQ mean = 1.02). Furthermore, 
CD18 and CD3e expression decreased dramatically (6 to 
10-folds reduction) in captured cells as compared with the 
whole mammary tissue, thus suggesting a negligible pres- 
ence of immune cells in laser-captured cells, further dem- 
onstrating the efficiency of LCM (Figure 4B). However, a 
weak expression was systematically recorded for Krtl4, in 
microdissected MEC, reflecting a slight contamination by 
mammary myoepithelial cells (MMC) during the course 
of the capture process, since it is known that Krtl4 is 
expressed by MMC only [20]. This is probably due to the 
very close proximity between MMC and MEC and to the 
difficulty of disassociating luminal MEC from MMC bor- 
dering the basal lamina that separates the epithelial layer 
from the extracellular matrix. 

MEC response during the first 24 h following IMI 
challenge with S. aureus 

Gene expression profiling of MFG at the early steps (be- 
tween 0 and 24 h) of the experimental IMI challenge with 
S, aureus was examined for the 5 goats of our experimen- 
tal design, using microarray technology to identify a reper- 
toire of differentially expressed genes in goat MEC in 
response to infection. Microarray results showed that the 
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Figure 4 Selectivity of MEC capture by Microdissection from infected mammary tissue, assessed by real-time qPCR. A) Section of 
infected mammary tissue witli tlie zone selected for laser capture microdissection (pink). B) Relative expression (± SEM) of specific marker for 
MEC iCSN3, green bar) and for putative contaminating myoepithelial cells (KrtM, red bar), macrophages and neutrophils (CD/8, orange bar) and 
lymphocytes {CD3e, yellow bar) in infected mammary gland (left) and in microdissected infected MEC (right), at 30 hpi. *Significance is relative to 
fold change in expression of infected MG (adjusted p-value < 0.05). 



most highly expressed genes, at the 4 time points of infec- 
tion (including 24 hpi), were the signature of lactating 
MEC gene expression. Indeed, genes coding for milk pro- 
teins {LALBA, CSNs and BLG) and ribosomal proteins in- 
volved in protein synthesis, RPL and RPS, gave strong 
signals. Likewise, amongst these highly expressed genes, 
we also found genes involved in synthesis and secretion 
of \fatty acids and triglycerides {FASN, XDH, ADRP, 
SCD and DGATl), as well as the gene encoding UDP- 
galactosyltransferase {B4GALT1) which participates in 
glycoprotein post-translational modification and lactose 
synthesis. Furthermore, microarray data are consistent 
with the results of qPCR analyses performed to assess 
the level of contamination of MFG by immune cells. In- 
deed, specific macrophage and neutrophil gene markers 
present on the sheep microarray indicate an absence of 
immune cells up to 18 hpi whereas the contamination 
by neutrophils and macrophages began to be noticeable 
at 24 hpi. Consequently, only the first three conditions 
(TO, T12 and T18) were considered in the differential 
analyses done to highlight MEC immune response to 
IMI. Globally, repeated measures ANOVA analysis of 
microarray data revealed that a Umited number of genes 
were impacted in their expression at the early stages of 
infection. Among the 15 208 probe set present on the 
Sheep Gene Expression Microarray (8 x 15 K), a total of 
47 probes (39 annotated genes) were significantly differ- 
entially expressed in response to IMI challenge after a 
Benjamini-Hochberg correction [28] to control the false 
discovery rate (adjusted ^-value inferior to 0.05). These 
genes (Table 2) showed at least a 1.5-fold difference in ex- 
pression in one of the pair wise comparisons (0 vs. 12 h, 12 
vs. 18 h). Among these 39 differentially expressed genes. 



only 12 genes are differentially expressed at 12 hpi, 
while all 39 are differentially expressed between 12 and 
18 hpi. This limited but increasing number of differen- 
tially expressed genes was consistent with the fact that 
gene expression analyses were performed at the very early 
stages of infection. Regarding the 0 h vs. 12 h comparison, 
fold changes were relatively low (ranging between 1.5 and 
2.7) with a majority of genes (11) down- regulated. Only 3 
genes were up-regulated corresponding to a chemokine 
(CCL-2), the immediate early response 3 protein (IER3) 
and the phosphodiesterase 4B cAMP-specific (PDE4B). 
Between 12 and 18 hpi, there were more differentially 
expressed genes, and the balance between up and down 
regulated genes was reversed. Only 3 genes were down- 
regulated during this period with low fold changes (1.6 to 
1.7). In contrast, 36 genes were up-regulated, with fold 
changes ranging between 1.5 and 21 (Table 2). Genes 
whose expression was highly increased, were mainly 
involved in the inflammatory response [interleukin-1 
receptor {IL-IRN), SlOO calcium-binding protein A12 
{S100A12), Interleukin-1 13 {11-1^)1 in chemotaxis [inter- 
leukin 8 {IL-8), chemokine C-C motif ligand 4 {CCL4)], in 
cell adhesion [tumor necrosis factor alpha-induced protein 
6 {TNFAIP6), platelet/ endothelial cell adhesion molecule 
(PECAMl)], in opsonization by enhancing phagocytic ac- 
tivity of infiltrating immune cells [pentraxin 3 {PTX3)] 
and in apoptosis [serglycin {SRGN), Bcl-2 related protein 
Al {BCL2Al)l 

Gene network analysis of differentially-expressed gene lists 

Systemic identification and grouping of differentially- 
expressed genes into biological networks was performed 
using the software packages Ingenuity Pathway Analysis. 
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Table 2 The 39 differentially expressed genes (adjusted p-value < 0.05 and fold change > 1.5) in Milk Fat Globules of 
goats (n = 5), in response to an experimental IMI with 5. aureus, in at least one of the two pairwise comparisons 
(TO vs. 12 hpi; 12 hpi vs. 18 hpi) 



Description 



Gene Biological function 
symbol 



Fold change T0/T1 2 Fold change T1 2/T1 8 
(sens of regulation) (sens of regulation) 



Arachidonate 5-lipoxygenase-activating 
protein 

Bcl-2-related protein Al 

Chemokine (C-C motif) ligand 2 

Chemokine (C-C motif) ligand 4 

C-type lectin domain family 4 
member E 

Catliepsin Z 

Early growth response protein 1 

Early growth response 3 

Fc fragment of IgG, low affinity Ilia, 
receptor 

Proto-oncogene protein c-fos 

Growth arrest and DNA-damage-inducible 
protein GADD45 alpha 

Immediate early response 3 

lnterleukin-18-binding protein Precursor 

Interleukin 1 beta 

lnterleukin-1 receptor 

lnterleukin-2 receptor gamma chain 

Interleukin 8 

Metallothionein-1 A 

Phosphodiesterase 4B, cAMP-specific 

Platelet/endothelial cell adhesion 
molecule 

Pentraxin-related protein 

Regulator of G-protein signaling 2 
Ras-related associated with diabetes 
S100 calcium-binding protein A12 

Protein strawberry notch homolog 2 

Serine dehydratase 

Alpha-1 -antitrypsin 

Solute carrier family 2, facilitated 
glucose transporter member 3 

Suppressor of cytokine signaling 3 

Spleen focus forming virus (SFFV) 
proviral integration oncogene 

Serglycin 

Tumor necrosis factor alpha 

Tumor necrosis factor alpha-induced 
protein 6 

Tristetraprolin 

Thioredoxin-interacting protein 



AL0X5AP Leukotriene metabolic process 

BCL2A1 Apoptosis 

CCL2 Inflammatory response, chemotaxis 

CCL4 Inflammatory response, chemotaxis 

CLEC4E Immune response 

CTSZ Proteolysis 

EGRl Regulation of transcription, cell proliferation 

EGR3 Apoptosis, chemotaxis 

FCGR3A Immune response 

FOS Regulation of transcription, cell proliferation 

GADD45a Cell cycle 

IER3 Inflammatory response 

IL18BP Regulation of transcription, cell proliferation 

IL-ip Inflammatory response, cytokine 
biosynthesis 

ILIRN Inflammatory response 

IL2RG Regulation of gene expression, alpha-beta 
Regulatory T cell differentiation 

IL-8 Chemotaxis 

MT-IA Cellular response to zinc and cadmium ion 

PDE4B Signal transduction 

PECAMl Cell adhesion 

PTX3 Regulation of phagocytosis 

RGS2 Cell cycle 

RRAD Protein transport 

S100A12 Inflammatory response 

SBN02 Regulation of transcription, cell proliferation 

SDS Cellular amino acid metabolic process 

SERPINAl Inflammatory response, acute phase response 

SLC2A3 Carbohydrate transmembrane transport 

S0CS3 Inflammatory response 

SPn Regulation of transcription, 
cell proliferation 

SRGN Apoptosis 

TNFa Inflammatory response, anti aptoptosis 

TNFAIP6 Cell adhesion 

TTP Inflammatory response 

TXNIP Response to oxidative stress 



/./ 

1.0 

1.6 (up) 
7.7 
1.3 

1.4 
1.1 

1.1 

1 .7 (down) 
/./ 

1 .8 (down) 

1.5 (up) 
1.3 

1 .5 (down) 

1.0 

1.2 

1.3 

1.4 

2.1 (up) 
1.4 

1.2 

1.2 

1.0 

1 .9 (down) 

1.4 
1.3 
1.3 
1.4 

1.4 

2.7 (down) 

1.3 

1.3 
1.0 

1.0 

2.1 (down) 



1.7 (up) 

5.1 (up) 

1.6 (up) 

6.8 (up) 

6.1 (up) 

1.5 (up) 

1.6 (up) 

2.3 (up) 
1.8 (up) 

2.5 (up) 

3.2 (up) 

2.7 (up) 

1 .7 (down) 

5.6 (up) 

9.8 (up) 

1.5 (up) 

8.5 (up) 

1 .7 (down) 

2.1 (up) 
3.0 (up) 

4.8 (up) 

2.2 (up) 

1.6 (up) 

8.0 (up) 

1 .6 (down) 

5.7 (up) 

2.1 (up) 

6.9 (up) 

2.4 (up) 
4.4 (up) 

8.4 (up) 

2.9 (up) 
9.1 (up) 

1.6 (up) 

1.5 (up) 
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Table 2 The 39 differentially expressed genes (adjusted p-value < 0.05 and fold change > 1.5) in Milk Fat Globules of 
goats (n = 5), in response to an experimental IMI with 5. aureus, in at least one of the two pairwise comparisons 
(TO vs. 12 hpi; 12 hpi vs. 18 hpi) (Continued) 



021030OOCX008067HT OOCX Ovisaries 
cDNA 


1,7 (up) 


4.0 (up) 


MGCl 65862 protein 


1.6 (down) 


3.4 (up) 


020605OCS41 1 01 6083HT 0CS4 
Ovisaries cDNA 


1.2 


21.0 (up) 


Taste receptor type 2 - Bos taurus 
(Bovine), partial (26%) [TCI 7737] 


1 .5 (down) 


2.7 (up) 



Adjusted p-value < 0.05 and fold change > 1.5 are in bold, in at least one of the two pairwise comparisons. 



Analysis of gene interactions and pathways was performed 
using the differentially-expressed gene sets (0-12 h and 
12-18 h) after staphylococcal stimulation, independently. 
We observed that top networks obtained at 12 hpi deal 
with cellular movement with 12 molecules on 14 mole- 
cules annotated taking part in this network. At 18 hpi, the 
most relevant biological networks were related to inflam- 
matory response and infectious disease with 27 molecules 
on 32 molecules annotated (Figure 5). Moreover, for each 
comparison, inflammatory response and cell to ceU sig- 
nalling and interaction were among the top biological flinc- 
tions that came up with highly signiflcant values (1.09 x 
10"^^-8.77 X 10"°^ and 1.09 x 10"^^-9.69 x 10"^^ at 12 hpi, 
respectively and 1.55 x 10"^^-8.04 x 10"^^ and 1.65 x 10"^^- 
7.43 X 10"°^ at 18 hpi, respectively. 

Using the Empirical Bayes Dynamic Bayesian Network 
(EBDBN) algorithm in order to identify temporal rela- 
tionships among the 39 differentially expressed genes, 
we were able to construct an interaction gene network. 
Using this approach, it appears that GADD4Say IL18BP 
and S0CS3 are highly regulated genes in the network 
with a high degree of outward connectivity to other 
genes that are for the most up-regulated (Figure 6). 
TNFa, PDE4B and S100A12 are conversely key genes en- 
coding molecules that activate or repress the expression 
of other genes of the network. 

Confirmation of differentially expressed genes by qPCR 

A total of 16 genes were chosen to be confirmed by 
qPCR and to substantiate the involvement of the key 
biological pathways identified. Together with the positive 
marker {CSN3) and contamination markers {CD18, CD68 
and GPR97), we chose eight up-regulated genes encoding: 
pentraxin 3 [PTX3), interleukins Ip and 8 [IL-ip, IL-8), 
chemokine 4 {CCL4), serglycin {SRGN), tumor necrosis 
factor alpha (TNFa), SERPINAl (SERPINAl), S100A12 
and four apparently non regulated genes encoding pro- 
teins known to be involved in the innate immune 
response: Toll-like receptor 2 {TLR2), cathelicidin 3 
(CATHLS), galectin3 {LGALS3) and the acute phase 
serum amyloid A3 {SAA3). All of these genes, with 
the exception of TLR2, that were found differentially 



expressed or not, as defined by the initial microarray 
screening, conserved the same status when analyzed by 
qPCR {P < 0.05). The reasons why TLR2 was not identi- 
fied in the microarray analysis whereas it unambigu- 
ously showed a differential expression when measured 
by qPCR (Figure 7) are not clear. This might be due to 
a weak specificity of the microarray probe, rather than 
to the use of a heterologous system since we did not 
find any difference between goat and sheep TLR2 
mRNA nucleotide sequences that match perfectly. For 
the other target genes, our results show the same trend 
(Figure 7) and correlate with microarray data, although 
systematically and significantly higher in qPCR. For 
instance, the expression of the gene encoding the pro- 
inflammatory cytokine IL-8 increased sharply between 
0 and 18 hpi (mean fold change = 94.7, with an import- 
ant variability between goats) whereas the specific 
MEC marker was constant and the contamination by 
infiltrating immune cells remained negligible (Figure 8). 
Indeed, we also observed a sharp increase in the expres- 
sion of different innate immune genes, such as PTX3, 
IL-ip, CCL4 and S100A12, at 18 hpi. Regarding the gene 
encoding the SAA3 acute phase protein, its expression 
was unchanged until 18 hpi whereas the expression of 
TLR2 showed a slight increase (fold change = 2.0). 

Taken together, these results strongly support the 
notion that during the early stages post-inoculation 
(before 24 hpi), MEC initiate a prompt response to 
infection and a sharp induction of innate immune 
genes that are predominantly associated with the pro- 
inflammatory response. 

Genes associated with the innate immune response are 
still expressed in MEC at 30 hpi 

Given that the MEC response to infection cannot be 
analyzed using RNA extracted from MFG at 24 hpi and 
beyond, due to MFG contamination by immune cells 
after the burst of SCC in milk, mammary tissue samples 
collected at slaughtering on goats at 30 hpi were used 
to monitor MEC contribution using the same qPCR ap- 
proach on RNA extracted from microdissected MEC. 
LCM is considered the gold standard method to analyze 



Brenaut et al. Veterinary Research 2014, 45:16 
http://www.veterinaryresearch.0rg/content/45/l/l6 



Page 1 2 of 20 




\ r. ^ ' 

1,505 



Direct interaction 



Indirect interaction 



0 


Complex 




Cytokine/Growth Factor 




Enzyme 


0 


Group/Complex/Other 


<§ 


Peptidase 


CO 


Transcription Regulator 


•if 


Transmembrane Receptor 




Transporter 


n 


Unknown 




Relationship 




Relationship 


Colors: 


■ 


down regulated 


■ 


up regulated 



Figure 5 Ingenuity Pathway Analysis (IPA) of the 39 differentially expressed genes in MFG transcriptome at 18 hpi vs. before infection 
(TO). The IPA legend defining the symbols depicted in IPA networks is given in the inset. Direct interaction is in solid line whereas indirect 
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Figure 7 Heat map of up- and down-regulated genes related to innate defense, selected from MFG transcriptome analysis of the IMI 
time course challenge with S. aureus. Gene expression was assessed using qPCR; genes sliown in red are up-regulated and tliose sliown in 
green are down-regulated in MFG infected at 12 hpi and 18 hpi, relative to MFG from the uninfected half-udder. Data are expressed in logio 
ratios with respect to the reference genes {RPS24 and PPIA), at each time point and represent the five biological replicates at each time point. 
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Figure 8 Expression of /L-8 in MFG, at different time points of 
infection (12 hpi and 18 hpi) induced by an S. aureus challenge. 

Abundance of mRNA arising from the gene encoding the 
pro-inflammatory cytokine IL-8 (yellow bar) is expressed relative 
(± SEM) to transcripts from CSN3, a specific MEC marker (green 
bar) and from CD18, a specific marker for macrophages and 
neutrophils (red bar). ^Significance is relative to the fold change 
in expression taking MFG from milk of the right uninfected 
half-udder as reference (adjusted p-value < 0.05). 



single cell-type gene expression, since it allows an isola- 
tion, under morphological control, of cells of interest 
(here MEC) in their physiological environment. How- 
ever, given the limited amount of material extracted, we 
only measured the expression of a few genes found up- 
regulated in qPCR at 18 hpi (see above). Thirty hours 
after the IMI challenge, all the genes analyzed displayed 
a higher expression than in uninfected MEC (Table 3). 
The expression of IL-8 and PTX3, for example, increased 
until reaching a fold change (322 ± 6 and 124 ± 5, re- 
spectively) of the same order as that found within the 
mammary parenchyma (418 ± 9 and 158 ±5, respect- 
ively), strongly supporting the notion that IL-8 and 
PTX3 are synthesized by MEC. This conclusion also ap- 
plied to SAA3, SERPINAl, and TLR2, thus suggesting 
an important contribution of this cell-type in initiating 
the inflammatory response. On the contrary, several 
genes found up-regulated at 18 hpi using MFG RNA, 
such as IL-ip, CCL4, TNFa and S100A12, showed a 
lower expression at 30 hpi as compared with the mam- 
mary parenchyma (Table 3). 

This under-expression at 30 hpi could be due either to 
a specific degradation of the cognate transcripts during 
the microdissection process or most likely to an evolu- 
tion in the gene expression patterns. The decreased ex- 
pression recorded for IL-l^ and TNFa at 24 hpi with the 
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Table 3 Changes in expression of some innate immune 
genes in the mammary gland parenchyma (MG) and in 
microdissected MEC, relative to non-infected tissue and 
non-infected microdissected MEC at 30 hpi, respectively 



Target 
genes 


Fold change in 
MG at 30 hpi 


Fold change in microdissected 
MEC at 30 hpi 


SAA3 


2.71 ±0.75 


3.01 ±0.70 


SERPINAl 


10.40 ± 1.02 


9.56 ± 0.89 


TLR2 


3.68 ± 0.87 


2.97 ± 0.53 


PTX3 


158.07 ±5.20 


123.83 ±5.11 


IL-8 


41 7.87 ±8.78 


321.85 ±6.55 


TNFa 


29.57 ± 1.77 


16.77 ±1.21 


CCL4 


8.30 ± 1 .26 


3.1 7 ±0.36 


S100A12 


42.41 ±1.02 


14.92 ±0.34 




32.32 ± 2.41 


7.55 ±0.57 



microarray analysis was consistent with this hypoth- 
esis, at least for 2 of the 5 goats (Figure 9). On the 
contrary, dynamics of the expression observed for 
S100A12 (Figure 9) and CCL4 (not shown), between 0 
and 24 hpi, showed a pronounced tendency to increase, 
which seems consistent with the type of effectors and 
their functions. Since S100A12 is also expressed by 
neutrophils and macrophages, this result is very likely 
due to the presence of leukocytes infiltrating the mam- 
mary tissue (mainly neutrophils) in the infected MG 
samples taken at 30 hpi. Interestingly, it should be 
mentioned that the gene encoding CD 14 antigen, a cell 
surface protein known to play a pivotal role in me- 
diating recognition of bacterial cell wall components 
[6,12], followed the same dynamic profile as that of 
TNFa in all five goats. When infection is well estab- 
lished, cellular factors recruited from the blood stream, 
the expression by MEC of pro-inflammatory cytokines 
and chemokines declined in favour of the expression of 
several acute phase proteins. 

Discussion 

It has proved to be difficult to reliably analyse in vivo 
the time course of the pathogen-specific immune re- 
sponse in the udder [29]. Several groups have designed 
and standardized experimental conditions to reprodu- 
cibly infect lactating udders from cows [6,30,31] and 
ewes [32] to generate meaningful results describing mo- 
lecular mechanisms of host-pathogen interactions. How- 
ever, infection of animals suffers from among-individual 
variation, especially regarding dairy ruminant species. In 
addition, the udder is a complex organ made of several 
cell types contributing differentially to its immune re- 
sponse [33]. Therefore, to characterize the MEC re- 
sponse, relevant model cells such as primary bovine 
MEC have been proposed for the molecular analysis of 
immune defense mechanisms in the udder [12,34,35]. 



Their potential to express genes of immune defense in 
response to stimulation with heat-killed bacteria was 
recently characterized showing pbMEC is capable of 
secreting large amounts of mediators of inflammatory 
response, namely IL-8 and TNFa in a pathogen dose 
dependent fashion [35]. Comparing the inflammation 
dependent regulation of factors contributing to the com- 
plement system between the udder and MEC to those 
from fully inflamed udders after infection with live 
E. coli pathogens also underlines some regulatory limita- 
tions of the pbMEC cell model, which does not properly 
reflect the mastitis induced regulation of chemokines 
and the complement system in the udder [14]. 

The present study is, to our knowledge, the first report 
describing the in vivo transcriptional host response of 
MEC in dairy ruminants, early after an IMI challenge 
with S, aureus. We assessed the contribution of MEC to 
initiate a pro-inflammatory response and to participate 
in the recruitment of neutrophils and macrophages 
at the beginning of staphylococcal infection. Usually, 
in vitro models are used to get an inside view of MEC 
function, but it remains a simplified and artificial model. 
A surrogate and less reductionist approach than the 
in vitro ceUular and invasive biopsy methods is the use 
of "cytoplasmic crescents" enclosed in MFC, first pro- 
posed by Maningat et al. [36] in humans and recently 
validated in goats [16]. However, RNA from MFG, which 
in theory only contain RNA from MEC, could be contami- 
nated by RNA from infiltrating immune cells, whose 
population increases in milk during the infection process, 
and which can be trapped in the cream layer along MFG 
isolation. Our results show that up to 18 hpi, RNA 
extracted from MFG is uncontaminated by RNA from 
immune cells (Figure 3). The expression level of neutro- 
phil and macrophage specific markers was 17 times less 
than in a non infected mammary gland whose immune 
cell population is not more than 5% in the cow mammary 
tissue [37]. Conversely, at 24 hpi, data should be inter- 
preted with care, and we considered it wiser not to treat 
MFG as fully representative of MEC due to the level of 
contamination of MFG by immune cells which can then 
not be negligible. RNA extracted from MFG reflects the 
specific physiological metabolic process in functional epi- 
thelial cells within the mammary gland, only until 18 hpi. 
However, regarding some genes known to be specifically 
and highly expressed in MEC, such as genes encoding 
milk proteins, RNA extracted from MFG remains repre- 
sentative of this cell type. As shown by qPCR experiments, 
LCM allows obtaining MEC enriched material with low 
contamination by MMC. In such a manner we were able 
to complete until 30 hpi our study of MEC response in 
their physiological context, even though this contamin- 
ation was slightly higher than that previously reported by 
Bevilacqua et al. [20] on healthy tissue. This is probably 
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Figure 9 Individual time-courses expression of selected genes assessed with RNA from MFG using microarrays, during the course of 
the IMI challenge by S. aureus. Selected genes are the following: S100A12, li-lfi, CDH, NFkB, IRAK4 and MyD88. Gene expression reflecting 
abundance of mRNA is given as fluorescence intensity (normalized signal). 



due to tissue damage associated with inflammation leading 
to a loss of the secretory epithelium cohesion (Figure 4A). 
However, since MMC were not shown to be able to in- 
duce the production of chemoattractants for neutrophils 
[38], we can assume that this slight contamination of 
MEC by MMC has in fact little impact on changes in the 
expression profile induced by the infection and therefore 
does not interfere with the contribution of MEC in innate 
immune gene expression. 

In addition, throughout the time course of the chal- 
lenge we confirmed that mRNA molecules encoding 
proteins involved in the primary function of the udder, i.e. 
milk synthesis, were not affected by the experimental 
infection, thus demonstrating that there is no significant 
loss of function in the mammary tissue, at least during 
the first 24 h. Therefore, the experimental model mea- 
sured changes in immune gene responses that are likely 
to be appropriate to study the initial stages of infection. 
On the contrary, it was clearly demonstrated in cows 
that signals are transmitted to the surrounding tissue 



and to the neighbouring uninfected quarters, very early 
(within the first 24 h) after contact with E. coli [39] . 

Transcriptome analyses provided here, even though 
performed with an incomplete repertoire of gene probes, 
revealed that gene expression profiles clearly changed 
at 18 hpi with a low number of differentially regulated 
genes [n = 39), concomitantly with a 5 x 10^-fold in- 
crease of S. aureus bacterial count in milk from infected 
udders. If similar changes occurred in response to the 
bacterial infection in the neighbouring uninfected udder, 
as reported by Jensen et al. [40] in bovine, then such 
changes may have masked some of the differences be- 
tween the infected and uninfected udder and this could 
be a reason why the list of differentially expressed genes 
is relatively short. 

The rather late change in gene expression, only 18 hpi, 
could be explained by the delay to reach sufficient bacter- 
ial concentrations in milk. Sensing the presence of bacteria 
is an important component of innate immunity [5]. 
Indeed, kinetic studies of experimental clinical mastitis 
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induced by the intra-mammary infusion of S. aureus in 
ovines showed that the inflammatory response (leucocyte 
reaction) is not initiated until bacterial concentrations 
reach 4.10^ bacteria/mL of milk [41]. Differentially 
expressed genes, identified in the microarray screen- 
ing, of which 15 were confirmed by qPCR, indicated 
that MEC respond to an IMI challenge with S, aureus 
by substantial increase in abundance of mRNA coding 
for a range of pro-inflammatory cytokines, chemotactic 
proteins, antimicrobial factors and acute phase proteins. 
Our results clearly demonstrate that MEC, which form the 
first line of defense against invading pathogens, are actually 
able to induce an up-regulation of immune-associated 
genes involved in the inflammatory response, thus substan- 
tiating previous in vitro studies [9,10,12,13,42]. We note 
that, as pointed out by Gunther et al. [14] in a comparison 
of primary culture of bovine MEC and udder, such a MEC 
model, though useful, is generally not considered to reflect 
what actually happens in vivo. It seems, however, that 
MFC provides more accurate information to analyze the 
MEC response, at least during the very first steps of an 
IMI challenge. 

Furthermore, some of the genes found up-regulated in 
MEC play a key role in orchestrating the regulation of 
other immune effectors. For instance, the temporal regu- 
latory network, built from the expression data using the 
Empirical Bayes Dynamic Bayesian Network (EBDBN) 
algorithm, underlines the first role of MEC in the im- 
mune response. This role consists of the recruitment of 
neutrophils through the up-regulation of PDE4B and 
TNFa and in the differentiation and proliferation of 
lymphocyte T through the up-regulation of GADD45a 
[43]. More precisely, it was observed that the gene en- 
coding the phosphodiesterase PDE4B, which is expressed 
in the bovine mammary gland [44], plays a key role in 
immune cefl recruitment. Ariga et al. [45] showed that 
neutrophil recruitment to the lung is impaired in PDE4- 
deficient mice, by inactivating the second messenger, the 
cyclic adenosine monophosphate (cAMP), thus abrogat- 
ing its negative effects on neutrophfl recruitment and 
activation. Jin et al. [46] provide evidence that the gene 
encoding PDE4B, which is expressed in mouse periton- 
eal macrophages, is involved in the control of TLR 
signaling, significantly reducing TNFa expression. 

The multifunctional pro-inflammatory cytokine TNFa 
was shown to induce the expression and release of IL-8 by 
bovine MEC [47]. Our results which are consistent with 
this, show that genes encoding TNFa and IL-lp and, to a 
less extent, IL-6 which are typical pro -inflammatory me- 
diators produced in response to TLR stimulation, are up- 
regulated at 18 hpi. IL-8, a potent chemokine capable of 
initiating an acute inflammatory response, plays a primary 
role in the recruitment of neutrophils into the gland [48]. 
Thus, upon recognizing the pathogen, MEC are able to 



send out rapidly a strong signal to recruit cellular factors 
of immune defense (macrophages, neutrophils) from the 
blood stream into the infected mammary gland. Indeed, 
this change in gene expression was then foUowed, 24 hpi, 
by a massive arrival of ceUs in milk (Figure 2B), mainly 
represented by neutrophils [49]. 

Our results also suggest that the gene suppressor of 
cytokine signafling 3 {S0CS3) occupies a crucial position 
in the network arising from the temporal analysis. This 
gene, important for the mammary tissue homeostasis 
[50,51], encodes an intracellular inhibitor of cytokine 
signaling that acts in a classical negative feedback loop 
[52,53]. This is consistent with the fact that the up- 
regulation of different cytokine encoding genes such as 
IL-ip, TNFa and CCL4 by MEC is restricted at the very 
early stages of infection, before a decrease at 30 hpi. This 
short-lived up-regulation was also observed during in vitro 
studies [15]. Nevertheless, in comparison with in vitro 
studies, after chaUenging primary bovine MEC cultures 
(pbMEC) with heat-inactivated preparation of S. aureus, 
the induction of cytokine-encoding genes such as IL-l^ 
and IL-8 occurred later on in our study. Gunther et al. 
[35] also observed the same results with a stronger and 
earlier increase in mRNA abundance in pbMEC chal- 
lenged with E. coli than that found in the udder as 
caused by acute mastitis. This ability to respond more 
quickly and strongly could be due to the fact that MEC in 
culture does not properly reflect the mastitis induced- 
regulation of cytokines in the udder and rapidly lose 
their functional features. The immune activity of MEC 
is regulated and balanced within the udder by factors 
locally delivered from other cell types [14]. These obser- 
vations reinforce the interest of studying MEC biology 
in its physiological context, which can be done using 
RNA extracted from MFG. This non-invasive technique 
allowed analyzing the in vivo contribution of MEC dur- 
ing the very first steps of infection on the same individ- 
ual, thus allowing a study of time-course variations 
while accounting for individual variability. 

On the contrary, increased expression of genes encod- 
ing molecules with bactericidal functions, such as 
S100A12 and PTX3, is the second very clearly identified 
immune function of MEC. PTX3 together with S100A12 
that could assist defense of the mammary gland against 
chronic and subclinical infections, have also been reported 
to be up-regulated in milk somatic cells in response to 
S. aureus infection in goats [54] and the resulting proteins 
were shown to be present in bovine milk, helping to re- 
solve the mammary tissue infection as well as potentially 
contributing to the maturation of the newborn calf epithe- 
lium and establishment of the newborn gut microbial 
population [42] . Concomitantly with the early induction of 
cytokine and chemokine-encoding genes, we identified 
an up-regulation of two genes coding for acute phase 
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proteins (APP): pentraxin 3 (PTX3) and a-1 and- 
proteinase (SERPINAl) which were also found in 
bovine milk at the beginning of mastitis [39,42]. We 
demonstrate here from LCM experiments that MEC 
contribute significantly to the increases of PTX3 and 
SERPINAl transcripts, even after the burst of immune 
cells in milk. Therefore, as previously proposed for 
SERPINAl [55,56] which is a potent inhibitor of IL-8- 
induced hematopoietic stem cell mobilization [57], 
these two acute phase proteins could be useful markers 
in the early diagnosis of inflammation. We also ob- 
served a slight up-regulation of the gene encoding 
antimicrobial peptide S100A12 which is known to be 
induced by inflammatory cytokines such as IL-1|3 [58] 
and which demonstrates chemotactic activity, attract- 
ing circulating leucocytes in inflammation conditions 
[59]. This protein, which probably amplifies the inflam- 
matory response via a recruitment of neutrophils [42], 
was reported to occur during all but the early stage 
response [60]. This is in agreement with our results. 
Indeed, it appears that MEC do not greatly contribute to 
its expression in the early stages of infection. However, the 
time-course of its response showed a clear tendency to- 
wards up-regulation from 24 hpi. Conversely, when the in- 
fection is well established, MEC are able to specifically 
express the acute phase protein SAA3 which can act dir- 
ectly on pathogens [61]. The gene encoding SAA3 is 
mainly induced by IL1-|3 and TNFa [62]. Studies have 
shown a growing interest for this protein and propose 
SAA3 as a good marker for mastitis. Interestingly, a high 
level of SAA has been observed in milk, in clinical [63] as 
well as in sub-clinical [64] bovine mastitis. Furthermore, 
whatever the dose and the invading pathogen, in vitro [65] 
as well as in vivo [17] studies reported high levels of 
mRNA molecules encoding SAA3 in response to infection 
or experimental challenge with E. coli or S, aureus cell wall 
components. However, 30 hpi we observed a very weakly 
(2.7 to 3-fold) enhanced expression of SAA3, This discrep- 
ancy could be explained, either by the fact that MEC in 
culture do not properly reflect MEC in their physiological 
context and/or that the udder is a complex tissue in which 
the dynamics of inflammation are different from that ob- 
served with a cellular model or even because the reading 
windows are not comparable. In previous studies (unpub- 
lished results), in agreement with the results reported by 
Eckersall et al. [64] in cattle, we observed that increases in 
acute phase SAA in milk of experimentally infected ewes 
occurred within 12 h after S, aureus infusion reaching a 
peak concentration at 72 h. Thirty hpi can still be consid- 
ered as the onset of the IMI challenge with S, aureus, and 
consistent with our results, Eckersall et al. [64] who 
recorded at the early stage of infection (48 h) the expres- 
sion of SAA3 (mRNA), reported a slight increase (3.6-fold) 
in abundance of mRNA molecules. To substantiate our 



results, it is worth noting that most of these proteins 
were found in proteomic studies performed from E, coli, 
S, aureus [66] and S, uteris [67] mastitis milk whey. In 
addition, it should be underscored that the pattern of 
the response to IMI is basically similar and the qualita- 
tive differences observed between species most likely re- 
flect different basal conditions and the type of invading 
pathogen [68]. 

Finally, as previously demonstrated in an in vitro study 
[12], we confirmed in vivo that the expression of TLR2y 
which is a key component in immune recognition of 
gram-positive bacteria by host cells, recognizing a wide 
spectrum of microbial components [69], was slightly 
affected by the infection (2.9-fold up-regulated 30 hpi) 
and essentially expressed by MEC (80% contribution). 
Thus, it is likely that MEC contain a fully functional and 
constitutively active Toll-Like Receptor signalling path- 
way that is slightly induced by the bacterial challenge 
but immediately responsive. Therefore, MEC have an 
intrinsic role in innate immune surveillance of the mam- 
mary tissue. The primary function is recruiting immune 
cells and the recognition likely occurs via the MyD88- 
dependent TLR signaling pathway. Stimulation of this 
pathway is considered the main mechanism enhancing 
expression of pro- inflammatory cytokines and activation 
of the innate immune response [70]. TLR2 stimulation 
triggers intracellular signalling cascades leading to the 
activation of NFkB which, in turn, leads to the activation 
of several genes encoding pro-inflammatory mediators 
such as TNFa, IL-lp and IL-6. In addition, TLR2 stimu- 
lation engages the production and release by MEC of 
chemokines such as IL-8, a potent chemo-attractant and 
activating factor of neutrophils, and, to a less extent, of 
cytokines CCL2 and CCL4 which display chemotaxis 
activity for monocytes and macrophages. Interestingly, 
due to a considerable among-goat variation, a number of 
genes possibly involved in MEC response to IMI chal- 
lenge with S. aureus, has not emerged from the micro- 
array analysis. Therefore, we analyzed some relevant 
genes, known to play crucial roles in the immune function, 
by comparing their gene expression profiles between 
goats, during the course of the infection. Thus, MyD88, 
an adapter protein involved in the Toll-Like Receptor and 
IL-IR signaling pathway, with a sharp increase at 24 h 
seemed to follow the same dynamic trends (Figure 9) 
as S100A12, whereas NFkB and IRAK4 (essential in the 
activation of NFkB), showed a strong increase at 18 hpi, 
before rapidly declining at 24 hpi, at least for 3 of the 5 
goats analyzed in this study. 

These data collectively demonstrate that we can easily 
access MEC biology in its physiological context via MFG. 
This non- invasive technique has allowed assessing the 
contribution of MEC during the first steps of infection, 
in vivo. Our results suggest that the MEC response to IMI 
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challenge involves both genes likely to affect pathogen 
function as well as genes (e.g. cytokines) that alter the 
behaviour of other cell types. Consistent with others, the 
picture emerging from this study is that the mammary 
epithelium is not simply a mechanical barrier but rather a 
functionally complex tissue capable of responding effect- 
ively to the intrusion of pathogens by altering its own gene 
expression profile and possibly that of other cell types to 
favour recruitment of immune cells and to synthesize bac- 
tericidal molecules. We report here a novel strategy of 
sampling to monitor the dynamics of gene expression in 
MEC which allows going further into the understanding 
of MEC immune capacity. Furthermore, modulation of 
these roles could be of importance in determining the out- 
come of an infection. The chronologically induced synthe- 
sis of cytokines at the inflammation site is important for 
pathogen clearance, wound healing and return to normal 
conditions. This approach could allow a better under- 
standing of MEC functions in animals showing different 
levels of genetic predisposition to mastitis, as tested with 
an in vitro study on sheep [71]. From a diagnosis point of 
view, several candidates for an early detection of mastitis 
were found. 

Despite incomplete and imperfect annotation (the probe 
with the highest differential expression unfortunately re- 
mains unknown), which may explain why some effectors 
known to be induced in response to infection were not 
identified in this study (e.g. epithelial p-defensins such 
as LAP of which the expression is delayed [72] as well 
as other antibacterial molecules), the repertoire of gene 
probes used (Sheep Gene Expression Microarray, 8 x 15 K) 
has nevertheless confirmed the involvement of a number 
of master cytokines and chemokines (IL-8, IL-1|3). Regard- 
ing IL-1|3, which is first synthesized as biologically inactive 
pro-IL-lp, its processing into mature, biologically active 
pro-inflammatory cytokine supposes activation by caspase- 
1 which remains to be reported for MEC, before being 
released in the extracellular milieu following a non- 
classical secretory pathway [73]. One cannot take for 
granted that mRNA detection guaranties protein ex- 
pression. The application of an original algorithm to 
construct gene networks of temporal regulation re- 
vealed the involvement of several factors known to play 
an important role in the inflammatory response, but for 
which we did not suspect that the MEC could be pro- 
viders. In a recent review [74] and two proteomic stud- 
ies [75,76], increasing of low abundant proteins such as 
IL-8, CD14, SAA, S100A12 and PTX3, were reported in 
mastitis milk from ewes and cows experimentally in- 
fected by E. coliy S, aureus and S, uteris or challenged 
with LPS, thus corroborating our transcriptional find- 
ings. Finally, the non-invasive sampling method (RNA 
extraction from MFG) provided an opportunity to per- 
form a dynamic study of IMI, and this in spite of the 



significant individual variability observed which can in 
such a manner be partly bypassed. It also gives the 
opportunity to achieve a large-scale vaUdation of the 
results in a significant number of individuals, which we 
plan to do in the near future. 
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